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Abstract 

The threshold behaviour of the cross section a(e + e~ — > t + t~) is analysed, taking into account 
the known higher-order corrections. At present, this observable can be determined to next-to- 
next-to-leading order (NNLO) in a combined expansion in powers of a s and fermion velocities. 



1 Introduction 



The Tau-Charm Factory, a high-luminosity (~ 10 33 cm -2 s _1 ) e + e~ collider with a centre- 
of-mass energy near the t + t~ production threshold, has been proposed as a powerful tool 

to perform high-precision studies of the r lepton, charm hadrons and the charmonium system 
[Q, |4|. In recent years, this energy region has been only partially explored by the Chinese BEBC 
machine (~ 10 31 cm -2 s -1 ). The possibility to operate the Cornell CESR collider around the t + t~ 
threshold || has revived again the interest on Tau-Charm Factory physics Q. 

A precise understanding of the e + e~ —* t + t~ production cross section near threshold is clearly 
required. The accurate experimental analysis of this observable could allow to improve the present 
measurement jj]] of the r lepton mass. The cross section a(e + e~ — > r + r _ ) has already been 



analysed to 0(a ) in refs. [g, ||, 10], including a resummation of the leading Coulomb corrections. 



The recent development of non-relativistic effective field theories of QED (NRQED) and QCD 
(NRQCD) |ll| has allowed an extensive investigation of the threshold production of heavy flavours 
at e + e~ colliders. The threshold bb [12, 13, |14| and ti ]l5[ production cross sections have been 
computed to the next-to-next-to-leading order (NNLO) in a combined expansion in powers of a s 
and the fermion velocities. Making appropriate changes, those calculations can be easily applied 
to the study of t + t - production. 

In this paper we will compile and analyse the known higher-order corrections to the t + t~ 
production cross section. Although some 0(a 4 ) contributions have not been computed yet, the 
dominant NNLO corrections can be already incorporated to the numerical predictions. One can 
then achieve a theoretical precision better than 0.1%. 

The perturbative 0(a 3 ) and 0(a A ) contributions are discussed in section |2[ Section ||| contains 
the relevant non-relativistic corrections at low velocities, generating 0(a n /v m ) effects. The photon 
vacuum-polarization and the initial state radiation contributions are accounted for in sections || 
and [|, respectively. In section |(| electroweak corrections are shown to be negligible. The numerical 
results for the e + e~ — * t + t~ cross section and our final conclussions are given in section [7|. Some 
technical details and detailed formulae are relegated to the appendices. 



2 The Perturbative Calculation up to 0(a 4 ) 

A NNLO analysis of a QED quantity, following perturbation theory in the number of loops, 
implies that contributions up to C(a 4 ) should be taken into account. Let us review the terms 
contributing to the total cross section of r production in e + e~ annihilation up to this order. 

At lowest order in QED, the r leptons are produced by one-photon exchange in the s-channel, 
and the total cross section formula reads 

* B (e + e-^T+r-) = 2 -^v(3-v 2 ) , (1) 



where v = y/1 — AM 2 / s is the velocity of the final r leptons in the center-of-mass frame of the 
e + e~ pair and M = m T is the r mass, v is an adequate expansion parameter for observables 
evaluated at energies close to the production threshold, since its value goes to zero as we approach 
this point. This makes a B vanish in that limit, being the global factor v in (|l]) of kinematic origin. 
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The quantum numbers of the r + t~ pair are those of the photon, J PC = 1 , which corresponds 
to allowed r + r~ states 3 Si and 3 Z?i in spectroscopic notation 2S+1 Lj. 

Electromagnetic corrections of 0(a) to a B arise from the interference between the tree level 
result and the following 1-loop amplitudes: 

i) 0(a) corrections to the e + e~7 vertex, 

ii) O(a) corrections to the t + t~7 vertex, 
Hi) vacuum polarization, 

iv) box diagrams (2-photon production). 

The contributions from i) and ii) are usually expressed in terms of the Dirac and Pauli form 
factors at one loop |Iq] . The corrections to the photon propagator Hi) are divided into two pieces: 
the leptonic contribution (I = e,/i, r), which can be calculated perturbatively in QED, and the 
hadronic contribution, where QCD corrections make a perturbative estimate at low energies un- 
reliable. The hadronic vacuum polarization can be related to the total cross section of hadron 
production by means of a dispersion relation. Finally, the interference of the tree-level amplitude 
with box diagrams iv) does not contribute to the total cross section, by virtue of Furry 's Theorem. 

Besides the above virtual radiative corrections, the cross section of 0(a 3 ) corresponding to the 
process of real photon emission, e + e~ — * t + t~ 7 , must be added. The Bremsstrahlung photon 
can be emitted by the initial or final fermion lines, but there is no contribution to the total cross 
section from the interference between both sets of diagrams, again due to Furry's Theorem. We 
clearly see that there is no overlap between initial and final state radiative corrections for the total 
cross section up to 0(a 3 ). A compilation of analytical expressions for all the terms mentioned 
above is found in Ref. |J. 

Let us consider next 0(a 2 ) electromagnetic corrections to the Born cross section. They come 
from several sources: 

• Interferences between the one-loop diagrams mentioned previously. The total cross section 
contributions from interferences between i), ii) and Hi) with box diagrams are again zero. 
The first term involving two-photon r production comes from the square amplitude of the 
box diagrams. 

• Interferences between the Born term and the following two-loop amplitudes: the electron and 
the r vertex two- loop corrections, contained in the expressions of the electromagnetic form 
factors, 0{a) corrections to the vacuum polarization, and three-photon production diagrams, 
for which only the real part is needed. 

• The 0(a 4 ) Bremsstrahlung cross section, coming from the interference between tree-level and 
one- loop diagrams with one radiated photon, and from tree- level diagrams with two photons 
attached in any of the fermion lines, corresponding to the process e + e~ — > t + t~ 77. It is no 
longer true, at this order, that initial and final state real radiation could not interfere. 

Recall that the spectral density ImII om (s) built from the electromagnetic current of the r leptons 
collects all final-state interactions, including both virtual and real radiation, for single-photon 
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production, that is 

48vt 2 q; 2 

a(e + e~ -» 7* -> t+ r~) = — ^ ImII om (s) , (2) 

where the tilde on <r distinguishes it from the physical total cross section which includes all kind 
of corrections. Relation @ results from a direct application of the optical theorem, and is more 
commonly written as the ratio 

R em (s) = 1 = 12vrImn cm ( S ), (3) 

i.e., normalizing a to the point cross section a p t = 4 g" 2 . The ratio R em is well suited for studying 
the non-relativistic dynamics of the r + r~ pair, as it fully contains the final-state interaction. 
Therefore, the threshold behaviour of the total cross section will be ruled by the expansion of R em 
at low velocities. The perturbative QED expression of R em is given in Appendix A up to NNLO in 
the combined expansion in powers of a and v. 

As long as we do not care about multiple photon production of r leptons, neither consider 
interference between initial and final state radiation, it is possible to factorize the total cross 
section as an integration over the product of separate pieces including initial, intermediate and 
final state corrections: 

f S , 1 

a[s) = / F(s, 



1 + e 2 IL 



(w 



2 

a(w) dw . (4) 



The radiation function F(s, w) |17J describes initial state radiation, including virtual corrections, 
and yfs is the total energy in the center-of-mass frame. The integration emerges to account for the 
effective energy loss due to photon emission from the e + e~ pair. As previously mentioned, (Q) is 
an exact relation for the total cross section only up to 0(a 3 ), but it includes the largest corrections 
coming from the emission of an arbitrary number of initial photons, which can sizeably suppress 
the total cross section. The 0{a i ) contributions not included in this analysis are those coming 
from two- and three-photon production diagrams, for which no velocity enhancement is expected 
in the threshold region and so represent pure 0(a 2 ) corrections ~ 0.005%, and the interferences 
between 2-photon Bremsstrahlung diagrams overlapping initial and final state radiation. However, 
we shall argue in section || that Bremsstrahlung contributions start at NNNLO in the combined 
expansion in a and v, and so they are beyond the scope of our analysis. 



3 Non-Relativistic Corrections: NRQED 

We now focus on the behaviour of the total cross section in the region just above the production 
threshold, where the small velocity of the produced r leptons is another relevant parameter, besides 
a. The final-state t + t~ interactions are encoded in the electromagnetic form factors. Written in 



terms of v, their expressions at one and two loops [18] show the existence of O (— ) and 0(\) power 



like divergences in the limit v —> 0. This is a general result for any number of loops: diagrams with 
n uncrossed photons exchanged between the produced leptons generate singular terms proportional 
to (^)", known as Coulomb singularities, which lead to a breakdown of the QED perturbative series 
in a when v — * 0. Resummation of such terms is therefore mandatory, and it was done a long time 
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ago [ 19 1 , resulting in the well-known Sommerfeld factor 



1^,(0)1* = - °f (5) 
1 — exp(— ait/v) 

multiplying the Born cross section (|l]). This factor corresponds to the wave function at the origin, 
solution of the Schrodinger equation, of two conjugate charged particles of mass M interacting 
through a Coulomb potential for positive energies E = Mv 2 . The appearence of this factor in the 
cross section can be intuitively understood, since the Coulomb interaction modifies the scattered 
wave function of the lepton pair. The 1/v behaviour of this factor makes the cross section at 
threshold finite. 

We clearly see that a NNLO calculation of the cross section in the kinematic region where 
a ~ v has to account for all terms proportional to v (a/v) n x [1; a; v; a 2 ; av ; v 2 ] with n = 1, 2, . . . 
The leading divergences (i.e. (^) n , n > 1) can be treated by using well-known results from non- 
relativistic quantum mechanics, but a systematic way to calculate higher-order corrections in this 
regime, such as x [v,v 2 , ■ ■ •], seems to be far from obvious, at least from the point of view of 
covariant perturbation theory in the number of loops. An adequate description would come from 
a simplified theory which keeps the relevant physics at the scale Mv ~ Ma, characteristic of the 
Coulomb interaction, allowing for a clear and systematic identification of leading contributions. 

NRQED [|TT[ was designed precisely for this purpose. It is an effective field theory of QED at 
low energies, applicable to fermions in non-relativistic regimes, i.e. with typical momenta p/M ~ 
»Cl. Interactions contained in the NRQED Lagrangian flB.l| ) have a definite velocity counting 
but propagators and loop integrations can also generate powers of v. With appropriate counting 
rules at hand, one can prove that all interactions between the non-relativistic pair r + r~ can be 
described up to NNLO in terms of time-independent potentials |2C[| , derived from the low-energy 
Lagrangian. It can also be shown that the contributions to the total cross section from diagrams 
with real photons emitted from the produced heavy leptons begin at NNNLO []. 

The key observable to study threshold effects in t + t~ production is the 2-point function, LI em (s) 
calculated at NNLO. Its fully covariant expression is written as the time ordered correlator of two 
electromagnetic QED currents of the r lepton = f^r: 



RUq 2 ) = % Im[-i J d 4 xe^ x (0\T (>(x) j M (0)t) |0>] . (6) 



4tt 

7 

Inserting the effective low-energy expression for the QED current, eq. flB.2|) , in the last equation, 
one can arrive to the basic relation between the spectral density at NNLO and the non-relativistic 



Green's functions |22j : 



C°(? 2 ) = % l™{Ci G(0, 0;E)-^ G c (0, 0; E)) , (7) 

with C\ a short distance coefficient to be determined by matching full and effective theory results. 
The details of this derivation are found in Appendix B. 

The Green's function G obeys the Schrodinger equation corresponding to a two-body system 
interacting through potentials derived from £ NR q E d at NNLO, that means suppressed at most by 

lr This result can be explicitly seen by going to the well-known expression for cr(e + e~ — > t~^j) at tree level (see 
e.g. |21|); the leading term is oc av 2 , i.e. NNNLO compared to LO terms ~ (a/v) n ~ 0{1). 
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a 2 ,a/M or 1/Af 2 , as dictated by the counting rules. Such potentials have been calculated in the 
literature [p3|, [24], p5l , and in configuration space they read 



VJr) 



V BF (r) 



+ 



r { \, An 

a(/x s ) x 2 

47T 



2/3i ln(/tr) + ai 



/3 1 2 (4 1n 2 ( / ur) + ^-) +2(2/5 1 a 1 + /3 2 ) ln(/}r) + a 2 



M 2 1 J + 2M 2 r 



2M 2 



S 2 (Sr 
-r - 3- — r- 



a(/i 3 



7T 



M 2 



SV 3 )( } 



(8) 



V 2 + — r (r V)V 



f (2S 2 -3)^)(r) 



3a(/x a ) „ 
+ 2M 2 r 3 ' 



(9) 



(10) 



^Ki(r) 



4M 3 



(11) 



Here a(/j, s ) denotes the electromagnetic coupling constant renormalized in the MS scheme at the 
scale fJL s = Hsoft- The latter is the renormalization scale set for the 0(a) and 0(a 2 ) corrections to 
the Coulomb potential (||), as determined in [23] and |24}| , respectively. Note that these corrections 
involve ultraviolet divergent light fermion loops (rut <C M), which cannot be accurately described 
within NRQED. The scale jl is equal to /U so/t e 7E , with 7 E the Euler constant, and the rest of 
coefficients in (S) take the values 



01 



n e . 



20 

— no . 
9 



02 



-4n e , 



«2 



55 \ /20 

y - 16 Ca h(+ y"< 



(12) 
(13) 



The constants 0\ and /?2 are the one- and two-loop coefficients of the QED beta function in the 
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MS scheme defined as 



din a 



0(a) 



01-^+02 
47T 



a V 
4tt/ 



+ 



(14) 



<iln ^i 2 

The number of active lepton flavors ri£ would be equal to two for interacting r's. If quark loops are 
included we should substitute ni — > re/ = (re^ + AcX^Q 2 )) being Q q the electromagnetic charge of 
the quark q (with mass lower than M). 

The Breit-Fermi potential V BF (see e.g. |25]j ) has been written in terms of the total spin i? 
and angular momentum L of the lepton pair. At NNLO, the heavy leptons are only produced in 
triplet S-wave states, so we just need to consider the corresponding projection of the V BF potential, 
(i.e. make S 2 = 2 and L = in (g)). V An is a NNLO piece derived from the first contact term 



written in £ NRQED , eq. (BJ.), which reproduces the QED tree level s-channel diagram for the process 
£ + £~ — > £ + £~ . In QCD this diagram connects qq color-octet states, so this piece is not present 
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in recent papers devoted to threshold electromagnetic quark production, where qq pairs can only 



be produced in color-singlet states. Finally, the term (11) is the first relativistic correction to the 
kinetic energy. 

The Green's function at NNLO will therefore satisfy the Schrodinger equation f] 

( " W ~ 4& + Vc{r) + VBF{r) + VUr) " E ) G(r ' r '' E) = 6(3){r " r ' ] ' (15) 

A solution of eq. (^) must rely on numerical or perturbative techniques. In the QED case, a 
significant difference between both approaches is not expected, being a such a small parameter |]. 
Consequently we will follow the perturbative approach, using recent results by Hoang, Penin and 
others fl3| , |2^ , pGf| , who calculated the NLO and NNLO corrections to the Green's function ana- 
lytically, via Ray leigh- Schrodinger time-independent perturbation theory around the known LO 
Coulomb Green's function: 

G{x,y;E) = G c (x,y;E) + 5G{x,y;E), 
5G(x,y;E) = - J d 3 z G c (x, z; E) (H - H ) G c (z,y; E) +••• 
= - J d 3 zG c (x,z;E)(-^+V BF (z) + V An (z) + VP(z) + VW(z)\ G c (z,y;E) +••• 

= 6 Ki , BF G + 5 An G + 6^°G + 5 2 G + ^ NLO G + • • • . (16) 

Here H = -V 2 /M + V^°(r) is the pure Coulomb Hamiltonian. We refer the reader to Appendix C 
for complete expressions of G c and the different 5G's, as calculated in the literature, and for a full 
discussion about the regularization procedure. Let us just quote here that the Sommerfeld factor 
(H), which appears in the LO cross section can be easily recovered from the basic relation (0), if 
one reminds the spectral representation of the Green's function 

a(r y, E) ^^^ + i^^m, (17) 

^ E n -E - le J (2vr) J E k - E - te 

with *$> n (r) the bound state's wave functions {E n < 0), and ^k( r ) corresponding to eigenfunctions 
of H with Ek = k 2 /M > 0. The LO spectral density is proportional to the imaginary part of the 
Coulomb Green's function, and so, from (|l7|), proportional to |^' CiB (0)| 2 , i.e. to the solution at the 
origin of the Schrodinger equation with the LO Coulomb potential. 

Finally, the short distance coefficient C\ must be fixed. The "direct matching procedure" [p?| 1 
allows a straightforward determination of Ci, by comparing the NNLO non-relativistic expression 
(0) with the result (|A.1| ) for R em , calculated in full QED keeping terms up to 0(a 2 ) and NNLO in 
the velocity expansion. The short distance coefficient C\ is then expressed as a perturbative series 
in a(fi hard ) 



2 Note that the Green's function built from the NNLO potentials also resums higher order contributions, like those 
diagrams with the insertion of more than one NNLO potential term. 

3 Although for heavy quarks the numerical solution of the Schrodinger equation has been shown to have more 
stable NLO and NNLO corrections, we should note that higher-order terms not under control are being resummed, 
some of which are cutoff-dependent pq |. 
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Cl(M, fJ, har d, Ufa 



1 + 



u(fJ>har 



7T 



-,(1) 



+ 



/ a(Hhard) \ 2 (2), , 

[ I C 1 {l^hardj l^fac) T ■ • • ) 



(18) 



where we have anticipated that Cj does not depend on any scale. The renormalization point 
Hhard, chosen for a-jjg in the short distance coefficient, needs not be equal to that governing the 
perturbative expansions of the correlators, fi sof t, which only contains long-distance physics Q The 
result of the matching reads 1 22] 

J 1 ) _ _A 



with 



7T^ 



.(2) 



39 



7T 



4 M 2 



1 M 2 
- ln 



f^har 



6 y-fac 



^(^-C3)+^ln2-l 



18 





[4/11 \" 




' 11 " 


+ 




+ n f 


[97T 2 \ 



(19) 
(20) 



The factorization scale \i fac is introduced to separate long and short distance contributions in the 
process of regularization (see Appendix C for details) . 



4 Vacuum Polarization 



We now turn over intermediate state corrections in formula (|j). For a complete NNLO de- 
scription of a(e + e~ — > r + r~), two-loop corrections to the photon propagator should be included. 
Despite having calculated the final state observable R em in the MS scheme, we can exploit the fact 
that the piece e 2 / [1 + e 2 II em (s)] is a renormalization group invariant, and so evaluate these set of 
corrections in the on — shell scheme, where decoupling of heavy fermions is naturally implemented. 
The on — shell renormalized vacuum polarization function is defined as 



nr(? 2 ) = n cm ( g 2 )-n em (o). 



(21) 



The light lepton contributions to the vacuum polarization are the standard 1- and 2-loop per- 



turbative expressions [28]: 

e 2 n e , M ( g 2 

with 

n<V) - 



^)n«(g 2 ) + QV)(g 2 ) + 0(a 3 



E - 

^ 3 

i=e,fi 



- - In 



q 2 \ 6m 2 
mi' q 



+ 



rt 2 \q 2 ) 



El 



£_4 C3 -ln(-il)-12^1„(-il) +0 (^ 
6 v mi J a 2 V mi / V q A 



(22) 

(23) 
(24) 



4 Differences are relevant when NNLO corrections are considered. 
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where we have only retained the relevant terms in the limit m 2 <C q 2 (mg are the pole light-lepton 
masses) . For the r contribution in the threshold vicinity q 2 > 4M 2 , resummation of singular terms 
in the limit v — * is mandatory. Under the assumption a ~ v, it is clear that we need to know NLO 
contributions to H T (q 2 ), which means retaining uniquely G c and <5f LO G in (|l6|), but performing the 
direct matching not only for the imaginary part but also for the real part (up to 0(a)): 

e 2 U^°(q 2 ) = ^ d ( G c (0, 0; E) + <5? LO G(0, 0; E)\ +ah 1+ a 2 h 2 . (25) 



The one- loop coefficient C\ was already obtained in (|i~9|), and h\, hi are fixed by demanding equality 
between ReII T calculated in full QED and expression (p5|). We get 

, 8 

1 / 21 \ 11 3 , 1 , M 

4ir z \ 2 J 32 4 2 fi fac 

In the hadronic sector, a perturbative estimate of the vacuum polarization in terms of free 
quarks is unreliable since strong interactions at low energies become non-perturbative. An alter- 
native approach consists in relating the hadronic vacuum polarization with the total cross section 
a(e + e~ — > 7* — > had), by using unitarity and the analyticity of II had (s): 



n had (s) = - / dt 



ImII had (i) 



t(t — s — ie) 

= r ^ a(e+e ~ ~ kad) . (27) 

167r 3 a 2 Jiml t(t — s — ie) 

Usually, a(e + e~ — > 7* — ► had) is conveniently parameterized and the unknown parameters fitted 
to experimental measurements or else related to phenomenological constants. In this paper we will 
make use of a parameter- free formula for a(e + e~ — * 7* — > had) in the low-energy region, where 
the non-perturbative effects are more important, and the perturbative result for the high energy 
part. Below 1 GeV, the electromagnetic production of hadrons is dominated by the p resonance 
(J PC: = 1 ) and its decay to two charged pions. The photon mediated tt + tt~ production cross 
section at a center-of-mass energy yfs is written as 

<?(e + e~ -> vr+vr-) = —1 1 - — | F(s) \ 2 , (28) 

with F(s) being the pion electromagnetic form factor defined as 

(^-1^10)=^)^--^+)^. 

In the isospin limit, only the 1 = 1 part of the quark electromagnetic current = Q u u^f^u + 
Qd d'j^d survives. An analytic expression for the pion isovectorial form factor was obtained in 
Ref. [29] using Resonance Chiral Theory pH] and restrictions imposed by analyticity and unitarity. 
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The so-obtained F(s), which provides an excellent description of experimental data up to energies 
of the order of 1 GeV, reads: 



F(s) 



M 2 p 



Mj-s- iM p T p (s) 



exp< 



96^/2 



ReA{ml/s,ml/M 2 ) 



(29) 



where T p (s) is the off-shell width of the p meson |}l] 



967T/2 



M ? S Im 



and 



A(ml/s,ml/M i 



967T2/2 

In (ml/M i 



A(ml/s,m 2 jM 2 p : 



+ 



8mi 



+ at In 



0~TT + 1 

ov - 1 



(30) 
(31) 



Formula (25) will be integrated in ( p7| ) up to an upper bound s p 
region above s p , we use the perturbative results of Imn had : 



e 2 imn U)d)S ( s ) = Yl N cQ 2 ^ 

q=u,d,s 



1 + 



7T 



1 GeV 2 . For the integration 

(32) 



for light quarks, in the zero mass limit, and 
e 2 Imn Cjfc ( S )= 9{s - Am 2 q )N c Q 2 q ^ 



q=c,i 



( 2m 2 
(1 + - 

V s 



4m 2 . 



+ ^C F 

7T 



+ 



2 4 

rrig m„/5 



+ 

s s 



9 

^2 



18 In — 2- 

s 



(33) 



for the b and c quarks Q In both (|32^) and (33) the first QCD loop correction to the quark vacuum 
polarization has been added, with a s the strong coupling constant. This simplified description is 
good enough to achieve an accuracy better than 0.1% for the e + e~ — ► r + r~ cross section. 

As a test of our method to calculate the hadronic vacuum polarization, we have computed its 
contribution to the running of a at the scale y/s = Mz , and compared it with the results of recent 
analyses devoted to this subject |3^, 33] . In the on-shell scheme the evolution of the electromagnetic 
coupling constant due to hadron polarization is commonly defined as 



a(s) 



a 



1 - A had a(s) 



5 At the energy scales of r production the b quark has not been considered in the particle content of the effective 
theory, but we will include it when running a to s — M%. The contribution of the top quark to (^) starts at 
y/s ~ 350 GeV, so it is highly suppressed by the t 2 factor in the denominator. 
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with 

A had a(s) = -47raRe [II had (s) -II had (0)] . 

At the scale Mz we get A had a(Mf ) x 10 4 = 268, to be compared with the values A had a(Mf ) x 10 4 = 
280 ± 7 and A had o(M|) x 10 4 = 276.3 ± 1.6, obtained in |§ and || respectively. Our simple 
estimate only deviates by 4% and 3% respectively, from those analyses. Considering that n had 
modifies a(e + e~ — > t + t~) near threshold by roughly 1%, our result has a global uncertainty 
smaller than 0.1% for the total cross section. 

Let us just mention that the theoretical description of the vector form factor of the pion has 
been improved in a recent paper 34] using a model-independent parameterization which can fairly 
reproduce experimental data coming from e + e~ — ► it + it~ up to higher energies, yfs ^1.2 GeV. 
With such results, we would gain knowledge on the hadronic contribution to vacuum polarization, 
but its numerical effect on our final estimate would not be relevant, considering the important 
features of the hadronic spectrum we are leaving out by using naive QCD perturbation theory from 

~ 1.2 GeV upwards. 



5 Initial State Radiation 



In this section we collect the radiative corrections to single-photon annihilation of the initial 
e + e~ pair. These include both virtual and real photon radiation, all of which are needed at 0(a 2 ) 
in a formal NNLO analysis of <r(e + e~ — ► t + t~). However, for the emission of soft photons (i.e. 
photons whose energy do not exceed an experimental resolution AE <C y/s), it is a well-known 
feature that the expansion parameter is not a but (a/ir) log(s/rrtg) log(E/AE), which may be 
quite large, making necessary to retain all terms of the expansion with respect to it. It is possible 
to perform such resummation by using an approach based on the Structure Functions formalism 



[17|. In this technique, the effect of initial state radiation is accounted for by convoluting the cross 
section without initial radiative corrections with Structure Functions for electrons and positrons, 
in analogy with a Drell-Yan process in QCD. In the leading logarithmic approximation (i.e. when 
only terms containing a factor L = log(s/m 2 ) with each power of a are retained) this formalism 
allows to represent the cross section in the form (Q): 



a(s) 



dx F(x, s) 



l + e 2 n(s') 



~a{s>) , 



(34) 



with the "available" center-of-mass energy after Bremsstrahlung loss defined as s' = s(l — x), and 
the radiation function 



F(x,s) = f3x 



+ ^ 



P-i 



a fir 2 1\ 3^ (3 2 fl T 2 37 
vr ( T ~ 2 ) + 4^ ~ 2l ( 3 T 



4(2 - x) In 



1 (l + 3(l-x) 2 ) 



ln(l 



6 + x 



(35) 



— (L 

IT 



1) 



The total cross section, at the kinetic energy above threshold E = y/s — 2M, is evaluated by convo- 
luting the photon-mediated cross section of t + t~ production without initial radiative corrections 
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with a weight function F describing such radiation effects, from an energy E down to E' ~ E — AE. 
The function F(x, s) becomes larger as x — > 0, i.e. for E' < E, and it strongly decreases as the x 
variable grows. 

Besides the leading (a/ir) n L n terms, expression ( |35"| ) also includes all 0{a) terms exactly. The 
analysis of higher-order terms, not included in F(x, s), is showed in |l7fl to go beyond 0.1% accuracy 
for the interval of energies 0.2 GeV< y/s < 10 GeV. We shall use eq. ( j35[) to evaluate initial state 
corrections to the total cross section. 



6 Electroweak Corrections 



The small corrections arising from r production through a Z boson can be easily incorporated 
in our basic formula (|34|). Electroweak production of heavy quarks including threshold effects has 
already been studied in previous papers [35, |3q] . The trivial part comes from the vector couplings 
of the Z current, which just add a term proportional to R em (s) to the total cross section: 



a 



■V* Z* 



(s) 



(7 



Pt 



1 



V e V T + 



2 i 2 

v R + a R 



where V£ and ai are the neutral-current couplings of charged leptons, 



■Rem (^) 



(36) 



-1 +4 sin 2 9 W 
4 sin 9\y cos 0w ' 
-1 

4 sin 9\y cos 6w 



(37) 
(38) 



At the t + t~ threshold, electroweak corrections are at least suppressed by terms of 0(8m^./M'^) ~ 
3 • 10~ 3 with respect to photon mediated production. Due to the further suppression induced by 
the couplings v e = v T ~ 0.05, these electroweak corrections represent a contribution below 0.0008% 
to the total cross section, and therefore, they will not be considered for our purposes. 

The non-trivial part of the electroweak corrections comes from the axial couplings of the Z boson 
with the non-relativistic final state fermions. For such contributions one needs to expand the QCD 
axial-vector current in terms of proper NRQED currents and then to construct the corresponding 
non-relativistic correlator, which is already a NNLO contribution describing the r + r~ system in a 
P-wave triplet state |35, 36]. However it is suppressed by 0(16m^/M|), so fully negligible in our 
analysis. 



7 Final Results for a(e + e 



T T 



We now use formulas collected in previous sections to analyse the behaviour of a(e + e~ — > t + t~) 
at threshold energies. Some of the parameters appearing in the different pieces take the following 
values: 



The r mass, extracted from [37|, is m T = 1777.03 ± 0.30 MeV. 
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Figure 1: The spectral density R em at low velocities in both QED and NRQED. 



• The two-loop running of the electromagnetic coupling constant, defined in the MS scheme, is 
needed to evaluate a(/j, SO f t ) and a(fj, hard ), which show up in the non-relativistic correlator and 
in the short-distance coefficient, respectively. The 1- and 2-loop coefficients of the /3-function 
were already given in (|l5). The reference value for the QED running coupling has been chosen 
by the relation a-^ (m^) = a, with a = 1/137.036 the ordinary fine structure constant. 

• The first QCD perturbative correction to the vacuum polarization of free quarks is propor- 
tional to the strong coupling constant a s (see eqs. (^) and (|33"D). At the energy scale of r 
production, it is appropriate to choose m T as the normalization point for a s ; the corresponding 
value is a s {m 2 T ) = 0.35 ± 0.03 |§. 

• The dependence on the various renormalization scales /x so/t , [i hard and fj, fac is very small. 
The most pronounced one comes from variations on the scale (J, so f t governing the combined 
expansion in a and v of the NRQED correlators. The logarithms of this scale over Mv, which 
show up in the non-relativistic Green's functions, suggest taking /i so/t ~ Mv ~ Ma ~ 13MeV 
to minimize the size of the NLO and the NNLO corrections. In fact, in the range 10 MeV 
< fJ-soft < 100 MeV the sensitivity to changes in this scale is reduced, and we have the smallest 
NLO and NNLO corrections to R em , varying in the whole range by less than 0.15% and 0.08% 
respectively. The residual dependences on the other two scales are fully negligible. 

The need for performing resummations of the leading non-relativistic terms (a/v) n [v,va,v 2 , . . .] 
is evidenced in Figs. |l] and ^. The spectral density R em , calculated in both QED and NRQED, is 
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Figure 2: Relative sizes of corrections to R em {s) as calculated in (a) NRQED and (b) QED. 



displayed in Fig. |l] as a function of the r velocity. The QED tree-level result vanishes as v — > 0, due 
to the phase space velocity in formula (|l|), which is cancelled by the first v _1 term appearing in the 
O(a) correction, making the cross section at threshold finite. More singular terms near threshold, 
v~ 2 ,v~ 3 , . . . arising in higher-order corrections completely spoil the expected good convergence of 
the QED perturbative series in the limit v — * 0. This breakdown is clearly seen in the behaviour 
of the 0(a 2 ) correction to the QED spectral density in Fig. [l]. This is no longer the case for the 
effective theory perturbative series, whose convergence improves as we approach the threshold point, 
as shown in Fig. |2|a, and higher-order corrections reduce the perturbative uncertainty inherent to 
any series truncated at a finite order. In the whole energy range displayed in Fig. ^a, the differences 
between the NNLO, NLO and LO results are below 0.8%, which indicates that the LO result, i.e. 
the Sommerfeld factor, contains the relevant physics to describe the threshold region, although 
NLO and NNLO corrections would be needed for more accurate descriptions of the total cross 
section. 

We can safely assume that the NNLO result for the spectral density has a theoretical uncertainty 
below 0.1% for energies close enough to threshold. At larger energies, the subleading contributions 
gain importance and the convergence of the double series in a and v is poorer, due to the higher 
powers of the velocity which are not taken into account. This is the opposite behaviour to that of 
the usual perturbative QED expansion, Fig. Qb, where the series convergence improves as we move 
far away the threshold. 

Adding the intermediate and initial state corrections we have a complete description of the total 
cross section of t + t~ production, as shown in Fig. ||. Coulomb interaction between the produced 
t's, governed by the parameter a/v, becomes essential right within few MeV above the threshold, 
and his effects have to be taken into account to all orders in this parameter, making the total cross 
section finite in this region. Initial state radiation effectively reduces the available center-of-mass 
energy for r production, lowering in this way the total cross section. We can verify that this 
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i i u i i i i i i i 

3.552 3.554 3.556 3.558 3.56 3.562 

Vs [GeV] 

Figure 3: The total cross section a{e + e~ — * t + t~) at threshold: at tree level (solid line); plus 
NNLO corrections to final state interaction, eq. M) (dashed line); and also including radiative 
corrections from the initial e + e~ state and from vacuum polarization, eq. (dash-dotted curve). 

reduction remains at higher energies above threshold by examining Fig. ||. A maximum energy for 
the soft photons, AE = 60 MeV has been chosen to perform the integration (fP|), 

We should emphasize that NNLO corrections do not modify the predicted behaviour of the 
LO and NLO cross section as calculated in previous works but are essential to improve the 

accuracy of experimental fits with higher precision data and, even more important, to guarantee 
that the truncated perturbative series at NLO gets small corrections from higher-order terms. In 
this way, we have shown that the theoretical uncertainty of our analysis of a(e + e~ — * t + t~) is 
lower than 0.1%, being the main sources of error our estimates of the hadronic contribution to 
vacuum polarization and of the initial state radiation. The former could be easily improved using 
similar techniques to those applied to estimate a(Mz), but at the energy point yfs = 2m T , including 
fits to <r(e + e~ — > hadrons) data, and the latter, being detector dependent, should be accurately 
monitored and their effects correctly implemented in data analyses. Nevertheless, the statistical 
uncertainty of the most recent experiments is still much larger than the theoretical one due to low 
statistics, and we should wait for future machines to improve it. 
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3.54 3.56 3.58 3.6 3.62 3.64 

Vs [GeV] 

Figure 4: Initial radiation effects in the total cross section o~(e + e~ — > t + t~) up to energies around 
100 MeV above threshold (v ~ 0.2); the dashed line represents a(e + e~ — ► t + t~), which does not 
include radiative corrections from initial state, as defined in equation ^j. 
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Appendix A 



The expression of the two-loop spectral density as calculated in perturbative QED and up to 
NNLO in the velocity expansion (y = y/E/M) reads J18|: 



nNNLO 
- ft 21oop QED 



+ 



a(/j, h 



7T 



3vr 2 



7T 



6v + —v 2 +0(v 



+ ct 2 (ji h ) 



vr 2 3 
8v + 2 



2+n/ -hi 



1, 4£ 2 M 2 



6 H-hard 



5 
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+ 



49vr 2 3 1 iVf 2 \ 

In — lnu w + 0(v ^ 



192 +2 K - 2 ^"vr 2 



f^hard 



(A.l) 



The constant k has already been defined in eq. (|20|). The renormalization point in the MS 1 scheme 
has been chosen equal to fj> hard , and M denotes the pole mass. 



Appendix B 

The NRQED Lagrangian relevant for our analysis reads 



£nrqed = ^(E 2 -B 2 ) + 



D 2 D 4 

iDf + CO h CA 7T + 

2M 8M 3 



2M 8M 2 V ' 8M 2 K ' 



+ " 



die 2 



6M 4 



(V>W 2X *) { X T °2<r{-hD) 2 i,) + h.c. 



+ .... (B.l) 



The lepton and antilepton are described by the Pauli spinors ip and x> respectively. Antilepton 
bilinears and higher-order operators have been omitted. The first line in eq. ( |B.l| ) is related to 
the kinetic term of the QED Lagrangian, with the bilinear ip terms coming from the expansion 
of the lepton relativistic energy up to 0(1/M 3 ). Second line terms reproduce the electromagnetic 
couplings of the leptons with photons of energy lower than M. Four fermion operators displayed 
in latter lines reproduce production and annihilation of an £ + £~ pair in a S-wave singlet (efe) or 
triplet (di and d$) state. Additional interaction terms between photon fields should be introduced 
to simulate fermion loops. The short-distance coefficients must be determined following the 
matching procedure up to a certain order in a, in order to absorb infinities arising in calculations 
beyond tree level. 

Which interactions are to be kept for a given precision (in a and v ~ p/M) is dictated by 
counting rules. The presence of two dynamical scales in the theory, the fermions three-momentum 
~ Mv, and their kinetic energies ~ Mv 2 , makes the NRQED counting rules more involved than in 
most effective field theories. While the factors of a in a specific diagram can be read off from vertex 
coefficients, powers of v are also generated by internal propagators and loop integrations. There has 
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been a hard discussion during recent years on how to organize calculations within NRQED/NRQCD 
in a systematic expansion in v |3{J, especially in the context of dimensional regularization. The 
situation seems to be clarified with the new formulation proposed in Refs. pC|, 41 1. In a cutoff 



scheme power counting rules for the velocity had been previously derived by Labelle [20] using 
time ordered perturbation theory together with the Coulomb gauge to separate the "soft" photons 
(with energy E~ ~ Mv) from the "ultrasoft" ones (-EL ~ Mv 2 ). Although quite troublesome for 
calculations beyond NNLO in the velocity expansion, these rules give the order in v of diagrams 
containing only soft photons by simple dimensional analysis. Following these rules one proves that 
the latter diagrams are all we need to describe low-energy interaction between the pair of fermions 
up to NNLO. Moreover, soft photons have an energy independent propagator and therefore all 
interactions up to NNLO can be described in terms of potentials, being this a highly non-trivial 
result which cannot be derived in the context of full QED covariant perturbation theory []. 

The effective jt + t~ coupling seen by the non-relativistic leptons is given by the expansion of 
the QED current in terms of the operators of the low-energy theory: 



(x) 



QM 2 



(^(-iB) 2 x )(,) + 



(B.2) 



We have only quoted the terms which are needed at NNLO. The first piece is a dimension-three 
current while the second has dimension-five and it is already of NNLO, as dictated by counting rules 
[20] due to the presence of the 1/M 2 factor. Notice that both pieces have quantum numbers 3 Si. 
There is another dimension five current, describing 3 D\ t + t~ production which however would not 
contribute to the NNLO cross section because the correlator of the product of a 3 Si current and a 
3 Z?i one vanishes. The Wilson coefficients of the NRQED 3 Si current encode the effects of the hard 
modes which have been integrated out. The coefficient b\ needs to be known at order a 2 , while 
&2 = 1 at NNLO. Inserting expansion ( [B.2D into the correlation function (^) leads to the NRQED 
expression of the ratio R em at NNLO 



47T 



Im(Ci Ai{E) 



1 



6M 2 



Co 



ME) 



where 



Ai = -i J (fx e iqx ( |r (V o- X ■ X j v i> 



A 2 



d x e 



iqx 







X ] ° 



■ip + h.c. 







(B.3) 



(B.4) 



(B.5) 



The short distance coefficients read C\ = (&i) 2 and C 2 = 1. The correlators A\ and A% contain the 
non-relativistic interactions derived from the NRQED Lagrangian. Such interactions, at NNLO, are 
purely described by instantaneous potentials, similar to those used in familiar quantum mechanics. 
Therefore, once the lepton pair is created by the external current with relative momentum k and 
until it is annihilated, the four point function describing their evolution reduces to a Schrodinger 
Green function for a two-body system with kinetic energy E = yfs — 2M, see Fig. |5[ The exact 

6 In terms of diagrams this statement means that only ladder diagrams with Coulomb-like photons and contact 
interactions with vertex factors up to NNLO contribute. Crossed ladder graphs vanish for soft photons. 
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Figure 5: Graphical representation of the NRQED vector- current correlator diagrams: the lepton 
pair £ + £~ is created and annihilated by the coupling ip^crx in (B.4), and all the intermediate 
diagrams of the £ + £~ non-relativistic NNLO interaction are resummed in the Green's function 
G(E). 



relation for A\ reads 



Ai(E) 



Tr 



6 



d 3 k f d 3 k 



(2tt) 3 J (2tt) 

lim G(r,r';E) 

;r'— »0 



■aG(k,k';E)a 



(B.6) 



where we have used the identity Tr(er -a) = 3 Tr(I) = 6. One can check that eq. ( |B,6D gives the right 
proportionality factor between A\ and G just considering the free case. There is no extra factor 
coming from the different normalizations of the relativistic and the non-relativistic quantities. 
For the A2 correlator we have 



A 2 (E) 



Tr 



-6 



d 3 k f d 3 k 



(2tt) 3 

'2 



(k 2 + k' 2 )G c (k,k';E) 



(2vr) ; 

Vi + V 2 r ,)G c (r,r';E)\ r y^ . 



(B.7) 



As A2 is already of NNLO, only the Green's function for the Coulomb potential shall be considered. 
Relation (|B.7| ) can be further simplified by using the Schrodinger eq. (|l5|), retaining just the LO 
piece of V c (r). For the imaginary part, we have 



^lmG c (r,r';E) 



(E — Vc°(r) ) Im G c (r, r'; E) = [E+- ) Im G c (r, r' ; E) , (B 



In the limit r,r' — * 0, the term a/rImG c represents an ultraviolet divergence which must be 
regularized. Following the direct matching procedure |2^] to fix the value of the short distance 
coefficient C\ allows us to drop power-like divergences, such as a/rlmG c \ r ^o, which must cancel 
with similar ultraviolet divergences in C\ in the final expression for the total cross section. Therefore 
we can safely substitute \a1A2 by Y1EM G c (r, r'; E)\ r y^, in (|B.3|) to get the complete relation 
between the spectral density at NNLO and the non-relativistic Green's functions: 



V 



(,n Jm^CtG^O-E) -i§G c (0,0;£ 



M 2 



3M 



(B.9) 
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where we have expanded the relation q 2 = (2M + E) 2 to first order in E/2M, which is already a 
NNLO contribution. 



Appendix C 



The well-known Coulomb Green's function [42|, solution of the LO Hamiltonian, at the origin 

. Mv \ T ( _ , a(fJ-s)' 



reads (v = y/E/M) 

G r c (0,0;E) 



4tt 



i v 



In 



-J +7 



W{ 1 



2v 



(C.l) 



where ^f(z) = -^logT(z) and T(z) is the Euler Gamma function. The superscript 'r' stands for 
'renormalized', since in the short distance limit r, r' — > the Coulomb Green's function, and some 



of the 5G, have 1/r and log(r) divergent terms. Following the lines of previous papers [22 



power-like divergences are subtracted and ultraviolet log-terms are regularized by introducing a 
cutoff fif and hence subtracting the energy-independent part. However, the imaginary part of G c 
has no ultraviolet divergent terms, so they would not contribute to the total cross section. This 
is not longer the case for the corrections 5 KitBF G and 8 An G, and their (imaginary part) residual 
dependence on the \if -scale will be canceled with the scale dependence of the coefficient C±, which 
is determined using the "direct matching procedure" [^7) described at the end of section ||. We 
quote the result for 6 Ki:BF G {Cp — ► l,Tp — > 1 and Ca — ► for the U(l) group) p2| : 
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(C.2) 



The integration for the V An potential is trivial, and the resulting (renormalized) correction 5 An G 
reads 



6 An G(0,0;E) 



- 2 %T^(0,0,£) 2 



(C.3) 



and [ 
read: 



The 0(a) correction to the Coulomb potential, Vc l \r), must be iterated twice because it is a 
contribution. The corresponding corrections <5^ LO G and 5i NIj °G have been calculated in [|1| 
respectively. The details of their calculation can be found therein. Their final expressions 

a(fj, s 
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and 



^ ' lm=0 
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V 2v 
The constants Cq, Cj 1 are defined in terms of fi\ (0): 

Co 1 = ai + 2/3 l7 , 
C\ = 2ft. 

The iteration of the 0{a ) piece, (r), was also computed in fl3[| : 

5 2 G(0, 0; J5) = (^1) 2 { £ F 2 (to) ((to + 1) (C% + L(v)Cf + L 2 (fc)C 2 2 
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+2 F ( m ) ( C o + L{v)C\ + (L 2 (v) + K{m))Cl - (2 7 + *(m + 1)) (cf + 2L(u)C 

m=0 

+L(v)Cl + (- 7 i» + lL 2 (t,)) Cf + i\r(u)Cf} , 
with the functions I(m), J(m, n), K(m), N(v) defined as 

I(m) = (m + 1) ^* 2 (m + 2) - #'(m + 2) + y - -— ^yy^ 

-2(*(m + 1) + 7), 

n+ 1 / T , , 1 \ 

J(m, n) = 2 *i(m - n) h 27 

m — n \ n + 1 / 

+2 (*(m -n + 1) - V(m + 1)), 

m — n 

K(m) = 2(*(m + 1) + 7) 2 + *'(m + 1) - ^ 2 (m + 1) + 27 2 , 
N(v) = L + ?pj L{y) - 7 L 2 (v) + ^L 3 (v) , 

and the constants 

2 

Cl = (y + 4 7 2 ) + 0l + 2(/3 2 + 2/3 iai ) 7 + a 2 , 

C 1 2 = 2(/3 2 + 2/3 1 a 1 ) + 8/3 2 7 , 
C 2 = 4/3 2 . 

None of the above mentioned Coulomb 5G corrections have energy dependent ultraviolet terms 
on their imaginary part, so no matching is necessary for them. 
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